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We carry out quantum inelastic scattering calculations of collisions of Rb atoms with inverting 
NH3 and ND3 molecules in the energy range between and 100 cm -1 , which are important for 
experiments using velocity-controlled molecular beams to probe scattering resonances. We focus 
on molecules initially in the upper level of the ammonia inversion doublet for j = 1, k = 1, which 
is low-field-seeking and can be controlled in a Stark decelerator. We calculate the integral elastic 
and state-to-state inelastic cross sections in the coupled states approximation. We demonstrate 
the presence of both shape and Feshbach resonances in the elastic and inelastic cross sections at 
low collision energies and discuss their origin in terms of the bound states of Rb-ND3 complex. 
We also consider elastic and inelastic cross sections in the ultracold regime, using close-coupling 
calculations, in order to assess the viability of sympathetic cooling of ND3 by Rb. The inelastic 
cross section for relaxation to the lower level of the inversion doublet is smaller than expected for 
such a strongly coupled system, but is still likely to be too large to allow sympathetic cooling for 
ND3 in low-field-seeking states. However, there is a good prospect that sympathetic cooling will be 
possible for molecules in high-field-seeking states, even when the collision partner is a magnetically 
trapped atom in a low-field-seeking state. 

PACS numbers: 34.20.-b,34.50.Cx,37.10.Mn 



I. INTRODUCTION 

Over the last few years, several methods have been de- 
veloped to cool molecules to temperatures below 1 degree 
Kelvin. These include buffer-gas cooling in cryogenic he- 
lium P, Q , Stark deceleration with switched electric fields 
[H 3 , velocity filtering [j| , optical deceleration with laser 
fields Q , and crossed molecular beam scattering pj . The 
availability of cold molecules opens up a new field of low- 
energy collision studies, in a novel regime where collisions 
are dominated by long-range forces and resonances. 

Stark decelerators can be used to control beam velocity 
as well as to reduce it Pulsed molecular beams can 
be generated with much smaller velocity spreads than is 
possible with conventional supersonic sources. This offers 
the opportunity to study scattering resonances at much 
higher resolution than has been possible in the past. Gili- 
jamse et al. @ have carried out a proof-of-concept exper- 
iment in which a velocity-controlled beam of OH radicals 
collided with Xe atoms in a jet. In this case the energy 
resolution of the experiment was limited to 13 cm -1 by 
the velocity spread of the Xe atoms. However, exper- 
iments arc under way to collide two velocity-controlled 
beams, which will provide much higher resolution. 

An alternative approach is to collide a beam with a 
sample of trapped atoms or molecules that are already 
nearly at rest. Sawyer et al. [Io| have recently measured 
collision cross sections for He atoms and H2 molecules 
in pulsed beams colliding with magnetically trapped OH 
radicals, and achieved a resolution of 9 cm -1 . Experi- 
ments to investigate the collisions of a velocity-controlled 
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beam of ND3 with trapped ultracold Rb atoms are also 
under way, and again should be able to provide much 
better velocity resolution. 

Scattering resonances in molecular collisions have been 
studied for many yea rs II ll a nd can be very important 
in chemical reactions [12i . Il3| | . In simple systems where 
only a few partial waves contribute, they may produce 
sharp structures in cross sections as a function of collision 
energy [l4l. [l5j. However, in more complicated systems 
with dense energy level patterns, the resonant structures 
may get lost in the background. It is therefore important 
to explore whether well-defined resonant structures are 
expected for collisions of molecules that can be velocity- 
controlled (such as ND3 and OH) and atoms that can be 
laser-cooled (such as alkali metal atoms). 

In this paper we study collisions of Rb atoms with 
ND3 and NH3 for collision energies between and 100 
cm -1 . We carry out quantum-mechanical calculations of 
integral elastic and state-to-state inelastic cross sections 
using the coupled states approximation and observe nu- 
merous scattering resonances. To understand the nature 
of the scattering processes, we study the resonance struc- 
ture of the scattering cross sections for individual partial 
waves. The resonances can be explained in terms of the 
bound states of the Rb-NH3 complex. We calculate the 
pattern of bound states near the lowest dissociation limits 
of the complex as a function of the end-over-end angular 
momentum and identify the bound states responsible for 
the resonances. 

There is a further reason for being interested in Rb- 
NH3 collisions. Methods such as Stark deceleration can 
slow molecules to velocities of a few metres per second, 
and the resulting molecules can then be trapped at tem- 
peratures of 10 to 100 mK 0, QjJ. It has not yet proved 
possible to cool such molecules further, towards the tem- 
peratures and phase space densities at which they might 
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undergo condensation to form quantum gases. There is 
great interest in methods that might be used to achieve 
this, and one of the most promising is sympathetic cool- 
ing [ID, [l9||, in which molecules are cooled by thermal 
contact with a laser-cooled gas of atoms such as Rb. 
However, sympathetic cooling can work only if the colli- 
sions are predominantly elastic rather than inelastic: if 
the molecules are in excited states and undergo inelas- 
tic (deexcitation) collisions, the kinetic energy released is 
usually enough to eject both collision partners from the 
trap. We therefore also carry out calculations of elastic 
and inelastic cross sections at the temperatures relevant 
to sympathetic cooling. 



II. THEORY 

A. Rb— NH3 interaction potential 

In previous work we obtained the potential energy sur- 
face for NH3 interacting with Rb from highly correlated 
electronic structure calculations [2(|. Spin-restricted 
coupled-cluster calculations including single, double and 
perturbative triple excitations [RCCSD(T)] were carried 
out on a grid of points in the intermolecular distance R 
and intermolecular angles 8 and x- The angle x corre- 
sponds to rotation of NH3 about its C3 axis. The an- 
gles Oi were chosen to be the points for 20-point Gauss- 
Lobatto quadrature, which include 8 = and 180° (with 
8 = corresponding to approach toward the lone pair of 
NH3). Calculations were carried out on a grid of R val- 
ues from 3.5 to 25 clq and for two azimuthal angles x — 
and 60°. 

For collision calculations, we need to evaluate matrix 
elements of the interaction potential between angular mo- 
mentum basis functions. This is most easily achieved by 
expanding the potential in renormalized spherical har- 
monics Cxfj,(0, x). After taking account of the C% v sym- 
metry of NH3 , the expansion takes the form 

v(R,e, x ) = E E tit* — 

A=0,1,... ^=0,3,6,... 1 + < W> 

x KV(#,x) + (-im,-^,x)]. (i) 

The dominant terms in this expansion are those with 
(i = and 3, and the corresponding coefficients may be 
written 



vx^R) 
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iVcos0)V/ 1 (-M)dcos0; 

(2) 

where Pxfj,(cos8) is an associated Legendre polynomial 
and 

V o (R,0) = l[V(R,6,0) + V(R,6,60°)} 
V 3 (R,8) = ±[V(R,9,0)-V(R,8,60°)]. (3) 

Vo can be viewed as the interaction potential averaged 
over x^ while V3 describes the leading anisotropy of the 



potential with respect to rotation about the C3 axis of 
NH 3 . To obtain the expansion coefficients in Eq. ([2]) 
at a given value of R for Rb-NH3, the potential func- 
tions Vfj,(R,8i) at each grid point 8i are first evaluated 
by reproducing kernel Hilbert space (RKHS) interpola- 
tion [21I . [22l | and the angular integrations are then carried 
out by Gauss-Lobatto quadrature. 

For calculations on Rb-ND 3 we need to reexpand the 
potential in a coordinate system based on the center of 
mass of ND3 instead of NH3 . This is done by first gen- 
crating a set of points for Rb-ND3 at the same grid of 
distances and angles as for Rb-NH3 . The center of mass 
of ND3 is shifted by a distance S towards the D atoms 
from that of NH3 . Each point is obtained by as 



V Rb - m3 (R,9,x) = V(R',8', X ) 



where 



R' = R{l + t 2 + 2icos( 
cos 8' = (cos 8 + t)R/R', 



Nl/2 



(4) 



(5) 



and t = S/R. The potential coefficients for Rb-ND3 are 
then obtained by quadrature exactly as for Rb-NH3, but 
with the transformed set of points. 

The Vq and V3 interaction potentials for Rb-NH3 are 
shown in Fig. [TJ The potential is strongly anisotropic; 
there is a very deep minimum, nearly 1900 cm -1 deep, 
at the N side of NH3 . This is due to the interaction of the 
lone pair of NH3 with the singly occupied 5s orbital of the 
Rb atom. There is also a shallow secondary minimum, 
only about 100 cm -1 deep, on the H side of NH3, arising 
from dispersion interactions. Both minima are on the C3 
axis of NH3. A detailed discussion of the alkali- atom- 
NH3 interaction can be found in the Ref. [2(| ■ 

The basis set which we used in our previous work (aug- 
mented triple-zeta for N and H and a basis set of triple- 
zeta quality for Rb) gave depths for the minima at lin- 
ear geometries of 1862 and 110 cm -1 (for N-side and H- 
side configurations, respectively). To improve this in the 
present work, we performed additional calculations for 
linear geometries near the minima with a significantly 
larger basis set. With the basis set extended to aug- 
mented quadruple-zeta for N and H, and with additional 
h functions on the Rb atom (with exponents 0.45,0.167), 
and an additional midbond g function (with exponent 
0.4), we obtained a global minimum well depth of 1881 
cm -1 with negligible (less than 0.1%) change in the equi- 
librium distance. All the potential points were then 
scaled by the ratio of the well depths calculated with 
quadruple-zeta and triple-zeta basis sets. The depth of 
the secondary minimum on the surface obtained by such 
scaling is 111 cm -1 , which agrees well with the value 
obtained from ab initio calculations with the quadruple- 
zeta basis set. In the scattering calculations described 
below we have used the scaled triple-zeta potential. 

For very large Rb~NH 3 separations the expanded po- 
tential is dominated by the vq.o and i>2,o terms, which 
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FIG. 1: The interaction potential of Rb-NH 3 from CCSD(T) 
calculations: Vq(R,6) component (upper panel) and Vs(R, 9) 
component (lower panel). Contours are labeled in cm -1 . To 
aid visualization, V3 is plotted only in the energetically acces- 
sible region defined by Vb < 0. 



V(R, 8, x) that couples the molecular rotations strongly 
to the intermolecular distance. The molecular rotations 
are in turn coupled to the nuclear spins in the molecule by 
hyperfine terms. However, for an atom in an S state the 
only terms that couple the molecular degrees of freedom 
to the atomic electron and nuclear spins are (i) very weak 
magnetic dipole-dipole interactions between atomic and 
molecular spins and (ii) any dependence of the atomic 
hyperfine coupling on the intermolecular distance R. If 
these very small terms are neglected, the collision prob- 
lem may be treated as the scattering of the molecule from 
an unstructured atom. 

The Hamiltonian of a rigid symmetric-top molecule 
may be written 

H mol = Bf + (C - B)j 2 (7) 

where j 2 is the operator representing the angular momen- 
tum of the molecule and j z denotes its projection onto 
the molecule- fixed symmetry axis z, while C and B are 
the rotational constants for rotation about the symmetry 
axis and an axis perpendicular to it. The cigenfunctions 
of this Hamiltonian are labeled by three quantum num- 
bers: j, the total angular momentum of the rotor; to, 
the projection of j on the laboratory-frame Z-axis and 
k, the projection of j on the body-fixed molecular sym- 
metry axis z, 

j 2 \jkm) = j(j + l)\jkm) 

jz\jkm) = m\jkm) 

jz\jkm) = k\jkm). (8) 



vanish as R~ e . Since the inherent error of the super- 
molecular approach is quite significant at large R, we rep- 
resent the long-range behavior of «o,o and «2,o using C^o 
and Cq_2 coefficients obtained from perturbation theory 
[20| . At short range we switch over to the full expansion 
using a switching function {1 + exp[a(i? — with 
a = 0.5 Oq 1 and b = 26 an. 



B. Scattering calculations 

The Hamiltonian that describes collisions between an 
atom and a molecule may be written 



2/' 



d 2 
dR 2 



R 2 



H a t + Hmol + V{ 



(6) 



where [i is the reduced mass for the colliding pair, L 2 is 
the operator for rotation of the collision partners about 
one another, and H ai and -£f mo i are the Hamiltonians for 
the isolated atom and molecule respectively. In the spe- 
cial case where either the atom or the molecule is in a 
closed-shell singlet state, the main part of the interac- 
tion potential Vlntcr is a single potential energy surface 



The N-H stretching vibrations are at sufficiently high 
frequency that they can be neglected in the present work, 
but NH 3 also has a bending (umbrella) vibration with 
two minima that arc equivalent in the free molecule. 
Tunnelling between these two minima produces a low- 
frequency splitting that corresponds to inversion of the 
molecule through a planar geometry. To describe this we 
introduce an additional degree of freedom, the inversion 
coordinate h. The vibration-inversion functions |+) and 
|— ) (corresponding to the lower and upper component of 
the inversion doublet) may be written 



|±) ~ [f(h eq - ft) ± f(h cq + h)] 



(9) 



where f(x) is a function that is sharply peaked around 
x = 0, so that /(/i C q ± ft) is a bending function that 
is sharply peaked around one of the equilibrium values 
=Fft e q (corresponding to one of the individual "umbrella" 
states of NH 3 ). 

NH3 occurs in two different nuclear spin configurations, 
corresponding to A\ and E symmetries. In the first case 
(referred to as ortho-NH 3 ) the molecule can occupy ro- 
tational levels with k = 0, 3, 6, . . ., while in the latter (re- 
ferred to as para-NH 3 ) only A; = 1,2, 4, 5... arc allowed. 
For ortho-NH 3 in k = states, only one component of 
each inversion doublet is allowed by symmetry. ND 3 is 
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slightly different and occurs in three nuclear spin symme- 
tries, Ai, A2 and E. Molecules in nuclear spin states of 
A\ and Ai symmetries can occupy rotational levels with 
k = 0, 3, 6, etc., referred to as ortho-ND3, while those in 
the E state can occupy levels with k = 1, 2, 4, 5, etc., 
referred to as para-ND3. Ortho- and para-species do not 
readily interconvert in collisions. By contrast with NH3, 
however, ND3 molecules in k = states can exist in both 
upper and lower components of the inversion doublet. 

The rotational functions adapted for permutation- 
inversion symmetry may be written [23[ 

[\jkm)T(-iy\j-km)]\±) (10) 

for para-NH3 and 

[\jkm) ± -km)}\±) (11) 

for para-ND3. If f(h cq — h) and f(h eq + h) do not over- 
lap significantly, evaluating of the matrix elements of the 
potential between the symmetry-adapted basis functions 
gives expressions that are isomorphic to those for the case 
of rigid symmetric top with parity-adapted monomer ba- 
sis functions. It is therefore sufficient to carry out the 
scattering calculation exactly as for a rigid symmetric 
top, but with the monomer energies for the even and odd 
combinations of \ jkm) and \j —km) shifted up and down 
from the rigid-rotor values by half the inversion splitting. 

The approach of restricting the basis set of inversion 
functions to just the lowest pair of states was introduced 
by Green (23|, |H] and Davis and Boggs , but was later 
verified to be a good approximation by Van der Avoird 
and coworkers [2(| [27| , who carried out calculations on 
Ar-NH 3 collisions to compare with theory that took the 
h coordinate into account explicitly. 

The rotational constants and inversion parameters 
used for NH3 and ND3 in the present paper are given 
in Table |U Since inelastic transitions between inversion 
states are of key interest for the present work, we allow 
the inversion splitting to be a function of rotational state 

a 

V = lS -Va[j{3 + l)-k 2 }-lS b k 2 (12) 

The energy levels for j = 1 to 3 for both ortho- and 
para-NH3 and ND3 are given in Table [II] In each case 
the zero of energy is the energy of the lowest allowed 
tunneling component for j = 0, k = 0. 

In the following sections we identify tunneling compo- 
nents by u (upper) or I (lower) instead of symmetry la- 
bels, because the u/l designation indicates the energetics 
more directly. 

The Rb-NH3 potential is very strongly anisotropic. 
Because of this, scattering calculations require large basis 
sets of rotational functions for convergence. In addition, 
we typically need to carry out calculations for many dif- 
ferent total angular momenta (partial waves) to obtain 
all the contributions to integral cross sections. We typi- 
cally need to include 100 to 200 partial waves at collision 
energies around 100 cm . 



TABLE I: Rotational constants [2Sl.[30| and inversion splitting 
parameters (with centrifugal distortions) [2f| used for NH3 
and ND3. Units are cm -1 . 



parameter NH3 ND3 

B 9.9441 5.1428 

C 6.2294 3.1142 

u 0.7934 5.337xl(T 2 

v a 5.05 xlO" 3 2.39 x 10" 4 

u b 1.998xl0" 3 9.61 x 10~ 5 



TABLE II: The energy levels of NH 3 and ND 3 molecules (up 
to j — 3) used in the present calculations. Units are cm" 1 . 
The labels o and p refer to ortho- and para- spin isomers. 
For NH3 in k — states, only the |( — inversion state is 
allowed for each j. 



state 


spin isomers 


NH 3 


l±> 


ND3 


l±> 


00 





0.0000 




0.0000 


+ 


00 









0.0534 




11 


P 


15.3871 


+ 


8.2570 


+ 


11 


P 


16.1735 




8.3100 




10 





19.1049 


+ 


10.2851 


+ 


10 









10.3385 




22 


P 


44.0305 


+ 


22.7424 


+ 


22 


P 


44.8058 




22.7949 




21 


P 


55.1837 


+ 


28.8282 


+ 


21 


P 


55.9499 




28.8802 




20 





59.6646 




30.8568 


+ 


20 









30.9102 




33 





85.1365 


+ 


43.4562 


+ 


33 





85.8969 




43.5080 




32 


P 


114.8786 


+ 


59.6850 


+ 


32 


P 


115.6145 




59.7356 




31 


P 


103.7254 


+ 


53.5992 


+ 


31 


P 


104.4704 




53.6503 




30 





118.5964 


+ 


61.7107 


+ 


30 









61.7641 





The most accurate approach for quantum inelastic 
scattering is to use close-coupling calculations, which ex- 
pand the total wavefunction for the collision system in 
a space-fixed basis set using a total angular momentum 
representation. In the present case the basis functions 
are labeled by quantum numbers j, fc, ±, which describe 
the monomer, and L, which describes the angular mo- 
mentum for rotation of the collision partners about one 
another. For each total angular momentum J and total 
parity, full close-coupling calculations include all possible 
values of L allowed by angular momentum coupling for 
each monomer level (j, fc, ±). However, such calculations 
give enormous basis sets for large J, and are prohibitively 
expensive except at the very lowest energies (where only 
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a few J values contribute to cross sections). 

A more affordable approach that is adequate at higher 
energies is provided by the coupled-states (CS) ap- 
proximation (3lL [3^ |. which was introduced for atom- 
symmetric top systems by Green [24|. In the CS ap- 
proach the scattering equations are written in a body- 
fixed frame and for each partial wave the basis functions 
are labeled by the monomer quantum numbers (J, k, ±) 
and the projection K (helicity) of j onto the body-fixed 
intermolccular axis R. If the centrifugal operator is ap- 
proximated by L(L+ l)/(2/ir 2 ). with L = J for all chan- 
nels, then the coupled equations factorize into indepen- 
dent sets labeled by L and K . The size of the resulting 
basis sets is independent of L, so that the total cost of a 
scattering calculation scales linearly with the number of 
partial waves included. The CS approximation is gener- 
ally expected to be accurate when the kinetic energy and 
the potential anisotropy arc both large compared to the 
rotational constant of the collision complex. 



C. Details of scattering calculations 

We carry out coupled-channel scattering calculations 
using the MOLSCAT package [H[ with the hybrid log- 
derivative/ Airy propagator of Manolopoulos and Alexan- 
der [2H . The inner starting point for the integration was 
chosen to be R = 4.2 ao, which is deep in the inner clas- 
sically forbidden region. Since the potential for 9 w 0° 
is very deep, the wavefunction oscillates very rapidly at 
short range. Thus in the inner region (i? < 13 ao) we 
use the fixed-interval log-derivative propagator [35| with 
a small step size of 0.015 ao- The log-derivative matrix is 
then propagated from 13 ao to 400 ao with the variable- 
step Airy propagator, which takes very long steps at long 
range and consumes only a small part of the total com- 
puter time. The Airy propagator was used with conver- 
gence cr iterion TOLHI = 10" 4 The resulting inelas- 
tic and elastic cross sections are converged with respect 
to all integration parameters to better than 1%. 

As mentioned above, the strong anisotropy of the po- 
tential energy surface causes slow convergence with re- 
spect to the basis set of monomer rotational functions, 
especially for ND3 , which has rotational constants about 
a factor of 2 smaller than NH3. In Table HTT1 we show the 
convergence of selected integral cross sections for Rb- 
ND 3 , together with relative computer times. The rota- 
tional basis set with j max = 21 and fc max = 8, containing 
210 monomer energy levels for para-ND3, gives a good 
compromise between accuracy and computational cost: 
typical calculations for a single energy and K value take 
~ 30 sec on a 2.0-GHz single-core Opteron processor, 
which means that calculations of cross sections converged 
with respect to the partial-wave expansion take approx- 
imately 2 hours per energy For Rb-NH 3 , a slightly 
smaller basis set with j max = 14 and fe max = 7 gave 
cross sections converged to ~ 1%. 

All calculations are for 87 Rb. 



TABLE III: The convergence of typical cross sections (in A 2 ) 
for RD-ND3 with respect to the maximum quantum numbers 
(jmax, fcmax) in the rotational basis set. The cross sections 
were calculated with the coupled-states method at an energy 
of 25 cm -1 with respect to j = 0, k = level. The computer 
time taken for a single partial wave, relative to the time for 
the (21,8) basis set, is also given. 



Basis 




0"llu->ll! 


f22«— ►llu 


time 


(15,5) 


1309 


62.22 


5.14 


0.2 


(18,5) 


1292 


59.42 


3.89 


0.3 


(21,5) 


1283 


62.67 


5.59 


0.5 


(15,8) 


1296 


57.39 


4.94 


0.3 


(18,8) 


1291 


63.11 


7.45 


0.6 


(21,8) 


1293 


59.70 


3.90 


1.0 


(15,11) 


1302 


61.83 


4.92 


0.5 


(18,11) 


1291 


62.09 


5.77 


0.9 


(21,11) 


1292 


59.42 


3.89 


1.7 



III. RESULTS 
A. Scattering cross sections 

The ND3 and NH3 molecules can be slowed by Stark 
deceleration in their low-field-seeking j = 1, k = 1 states, 
which correlate at low field with the upper level of the in- 
version doublet. Since ND3 has a much smaller tunneling 
splitting than NH 3 , its Stark effect is more nearly linear 
and it is easier to decelerate. We have calculated the en- 
ergy dependence of the state-to-state integral cross sec- 
tions for molecules initially in the upper inversion state 
of ND3 and NH3 for collision energies between and 100 
cm -1 . The results are shown in Fig.O 

For the whole energy range considered, the elastic cross 
sections an u ->iiu are large compared to the total inelas- 
tic cross sections for both ND3 and NH3 . Fig. [2] includes 
the elastic cross section obtained from the scmiclassi- 
cal background formula [3]} for a pure R~ 6 potential 
(dashed red line). For both Rb-NH 3 and Rb-ND 3 the 
elastic cross sections follow the background formula fairly 
closely, with slow glory oscillations superimposed on the 
background. The total cross sections (elastic + inelastic) 
are actually in even better agreement with the semiclas- 
sical model than the elastic cross sections. There is also 
structure due to scattering resonances at low collision en- 
ergies, but it is not very strong with respect to the back- 
ground for the elastic cross sections. At higher collision 
energies the resonances are lost in the background. 

The overall magnitude of the inelastic cross sections is 
at first sight less than expected. In a strongly coupled 
system where every collision that crosses the centrifugal 
barrier leads to inelasticity, the total inelastic cross sec- 
tion is given by the Langevin capture formula [38j . 

fc 6 \ 1/3 

^capture = 37T I — J . (13) 
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FIG. 2: Elastic and total inelastic cross sections for RD-NH3 
(upper panel) and RD-ND3 scattering (lower panel) from the 
upper component of the inversion doublet for j = 1, k = 1 
state. The smooth dashed lines show the result of the semi- 
classical background formula for elastic cross sections and of 
Langevin capture theory for the total inelastic cross sections. 



FIG. 3: Partial-wave contributions to inelastic Hit —* 111 
cross sections for Rb-ND3 at collision energies of 0.6 cm" 1 
(upper panel) and 10 cm -1 (lower panel). The straight dashed 
lines show the result corresponding to Langevin capture the- 
ory, with a cutoff determined by the heights of the centrifugal 
barriers calculated for the full potential. 



The Langevin result is shown as a dashed green line in 
Fig. [2] and it may be seen that the actual llu — > 11/ 
inelastic cross section is substantially below it. This is 
surprising in view of the very large anisotropy which di- 
rectly couples the 11/ and 11m states. Further insight is 
given by the partial- wave contributions to the llu — * 111 
cross sections, which are shown for collision energies of 
0.6 cm -1 and 10 cm" 1 in Fig. [31 In this case the capture 
cross section corresponds to an elastic S-matrix element 
of zero for every L below a cutoff due to the centrifugal 
potential. The partial inelastic cross section is far less 
than the capture value across most of the range of L. This 
contrasts with the behavior observed in other strongly 
coupled systems, such as high-energy HF-HF collisions 
39 1, rotationally inelastic collisions in Ar-N2 at 300 K 
40l | . and barrierless reactions in alkali metal atom + di- 
atom collisions at energies above 1 mK [41], [42|, EH, [44[ , 
where the partial cross sections are close to the capture 
value for all L. 

A possible explanation for the low inelasticity for Rb- 
ND 3 is that, despite the large anisotropy, the collisions 
are approximately adiabatic in the rotation-inversion co- 
ordinates. Fig. [?] shows "adiabatic bender" curves for 



Rb-ND3, which are eigenvalues of the rotation-inversion 
Hamiltonian for para-ND3 at fixed R. It may be seen 
that the curves correlating with the llu and 11/ states 
stay far apart for all values of R, with no avoided cross- 
ings between them where strong nonadiabatic transitions 
would be expected. 

The resonances are much stronger in the inelastic cross 
sections than the elastic cross sections. They are partic- 
ularly strong for collision energies up to about 20 cm" 1 , 
as can be seen in the expanded view of the energy de- 
pendence of criiu->iii hi Fig. [5] For Rb-ND3 at low en- 
ergies the resonances can enhance inelastic cross sections 
by up to a factor of 2. As the collision energy increases 
the resonances again become weaker and wash out. We 
see significantly more resonances for Rb-ND3 than for 
Rb-NH 3 . A detailed description of the origin of the res- 
onances for low kinetic energies will be given in section 
III.B. 

The coupled states approximation is not expected to 
be accurate at very low energies, because it approximates 
the centrifugal terms in the potential and neglects cou- 
plings between different hclicities. Unfortunately, full 
close-coupling calculations for higher energies are pro- 
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FIG. 4: Adiabatic bender curves for Rb interacting with para- 
ND3. The line types indicate the level that each curve corre- 
lates with at long range and not the character of the adiabatic 
state at short range. The cross indicates the position of the 
absolute minimum of the potential surface and the horizonal 
line indicates the energy of the lowest bound state, —1658.252 
cm -1 . 
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FIG. 5: Resonances in the inelastic 11m — > 111 cross sections 
for Rb-NH3 and Rb-ND3 at low collision energies. 



hibitively expensive for partial waves corresponding to 
large values of J. For example, the basis set used above 
(jmax = 21 and /c max = 8), which gives 210 channels 
for para-ND3 in the CS approximation, produces 2745 
channels for each parity in CC calculations for J > 21. 
However, it is possible to make the comparison for a re- 
stricted basis set, and the results are shown in Fig. [5] for 
a basis set with j max = 6, fc ma x = 5. The terms that 
are neglected in the CS approximation cause small shifts 
in resonance positions and intensities, so that it is not 
possible to make a valid comparison at any individual 
collision energy. However, the general extent of the reso- 
nance structure is similar in the two calculations and the 
inelastic cross sections typically agree within 5% above 
30 cm -1 . We can therefore have confidence in the general 
features of the CS calculations for the full basis set. 
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FIG. 6: Comparison of elastic Uu — > Hit (upper panel) and 
inelastic llu — > 11/ (lower panel) cross sections for Rb-ND3 
calculated with CS and CC methods for a restricted basis set 

With j max = 6, fcmax = 5. 

The state-to-state cross sections from the llu state to 
other j = 1 ... 3 states are shown as a function of colli- 
sion energy in Fig. [7] As the transition to each additional 
rotational state becomes energetically allowed, the corre- 
sponding excitation cross section rises sharply from zero. 
This behavior could be seen readily in experiments that 
use state-selective detection for the scattered ND3 . How- 



8 



ever, since the cross sections for rotational excitation are 
very small compared to those for 11m — > 111 relaxation 
at these energies, we do not see any significant changes 
in the total inelastic cross section as new channels open 
up. 
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FIG. 8: Near-threshold behavior of individual partial-wave 
contributions to the llu — ■» 21? inelastic excitation cross sec- 
tion, from close-coupling calculations, as a function of the 
excess energy above the threshold for excitation. 
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FIG. 7: State-to-state inelastic cross sections from the 11m 
state of Rb-NH3 and Rb-ND3 as a function of collision en- 
ergy, from coupled-states calculations. 

Fig. |8] shows the threshold behavior in more detail, 
for the contributions to the cross section <7n u ^2ii 

ob- 
tained from CC calculations for different partial waves 
J with the full basis set. CC calculations are feasible in 
this case because only low partial waves contribute just 
above threshold. At very low energy, each partial-wave 
contribution follows the Wigner threshold law [45[ and 
is proportional to (S con - £thresh) Ltinal+1/2 , where E coU 
and E'thresh are the collision energy and the threshold en- 
ergy respectively. L takes the lowest allowed value for 
the threshold level: for jfi na i = 2, L = for J = 2, 
L = 1 for J = 1 and 3 and L = 2 for J = and 4. For 
any final j there is always one total angular momentum 
(J = j) and parity combination that allows L = in the 
outgoing channel, and this dominates the integral cross 
section for excess energies below 10 -4 cm" 1 . Above this, 
however, additional partial waves contribute and simple 
(£coii — -E'thrcsh) 1 / 2 behavior is not expected for the inte- 
gral cross section. 



For Rb-NH.3, the llu — > 111 transition is much 
stronger than any other inelastic transition for the whole 
range of collision energies considered here. For Rb-ND3, 
ciiu^2iz becomes comparable to <J\i u ^iii for kinetic en- 
ergies around 100 cm . The resonances in individual 
state-to-state cross sections are in some cases quite strong 
compared to their backgrounds, especially closely above 
the thresholds. However, since their magnitude is much 
smaller than <7ii u ^iiz, these resonances are not visible in 
the total inelastic cross section. 

There are clear propensity rules that apply to the cross 
sections, as observed previously for He-NH.3 [23], 0, [H, 
HU and Ar-NH 3 [H,^ collisions. The general effect has 
been discussed by Alexander [48| . For fc-conserving colli- 
sions, inversion-changing transitions are preferred when 
Aj is odd and inversion-conserving transitions are pre- 
ferred when Aj is even. Thus <7n tl _»2ii is larger than 
ciiu^2iii, since the 11m and 21/ channels are directly 
coupled by V10 potential term. Conversely, for collisions 
with Ak = 1, which are driven by the V33 potential term, 
inversion-conserving transitions are preferred when Aj is 
odd, whereas inversion-changing transitions are preferred 
when Aj is even. Thus an u —>22u is much larger than 
o\\ u ^22i ■ It is interesting that these propensity rules 
survive in a system as strongly coupled as Rb-NH3. 



B. Analysis of scattering resonances 

There are in principle two types of resonance that 
might be seen in low-energy collisions. Shape resonances 
correspond to quasibound states that are confined be- 
hind a centrifugal barrier, whereas Feshbach resonances 
correspond to quasibound states that reside principally in 
channels corresponding to internally excited states. The 
height of the centrifugal barrier for each L is given ap- 
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proximately by [4^ 



V, 



L 

max 



h 2 L(L + 1 



, 3/2 



/' 



(54C 6 ) 



-1/2 



(14) 



This function is compared with the actual barrier heights 
obtained from the adiabatic bender curves for Rb-ND3 
in Fig. [9] it may be seen that the long-range formula 
(fT4| is accurate for low L but overestimates the bar- 
rier height by about 10% by L = 30. However, since 
Rb-ND 3 has a strongly attractive long-range potential 
{C&o = 523 Ehdo) and large reduced mass, the heights of 
the centrifugal barriers arc much smaller than in lighter 
systems. As will be seen below, most of the resonances 
observed in the present work are Feshbach resonances. 
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FIG. 9: Heights of the centrifugal barrier as a function of L 
predicted by long-range formula (|14|l and obtained from the 
full potential. 

The origin of the resonances can be understood if we 
study the contribution to the inelastic cross sections from 
individual terms in the partial-wave expansion. The top 
two panels of Fig. [TU] show the Uu — ► III partial cross 
sections for L = 0...6 and L = 7... 11 for Rb-ND3 
collision energies up to 2 cm^ 1 , where there are very 
strong resonances. The bottom panel shows the S-matrix 
cigenphasc sums (49| for L = 7 ... 1 1 , which show a sharp 
rise through it as the energy passes through a resonance. 

In general each partial inelastic cross section shows a 
non-resonant peak near the corresponding barrier max- 
imum and then dies off at higher energies. The non- 
resonant peak is quite sharp for low L but broadens as 
L increases. In particular, the peaks below E co n = 0.1 
cm -1 for L < 6 are non-resonant. Superimposed on the 
non-resonant background are peaks and troughs due to 
resonances. Below the barrier maximum, both shape and 
Feshbach resonances may occur; however peaks due to 
Feshbach resonances are suppressed when they are below 
barriers and it is usually shape resonances that give large 
peaks in this region. 

Above the barrier maximum for each L, Langcvin cap- 
ture theory would predict a partial inelastic cross section 
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FIG. 10: Contributions from partial waves L — ... 6 (top 
panel) and L = 7... 11 (center panel) to the inelastic llu — » 
111 cross section for Rb-ND3. The resonances are labeled by 
Roman numerals corresponding to the series of bound states 
shown in Fig. 1111 The bottom panel shows eigenphase sums 
for L = 7 ... 11. The dashed vertical lines in the bottom panel 
indicate the positions of centrifugal barriers. 



Except at resonances and occasionally near the barrier 
maximum, the partial cross sections are generally much 
smaller than this. It is thus evident that low-energy in- 
elastic scattering in Rb-ND3 is dominated by resonant 
effects. 

We also carried out bound-state calculations close to 
the j = 1, k = 1 thresholds to help identify the states 
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FIG. 11: Near-threshold bound states of Rb-ND 3 (full circles) 
for helicity K = 1 as a function of L(L + 1). Open circles in- 
dicate the positions of the strongest resonances in Fig. 1101 
The arrow points out an example of an avoided crossing be- 
tween bound states. The horizontal solid lines at 8.257 and 
8.310 cm -1 show the 111 and 11m thresholds, respectively. 
The dashed line shows the position of centrifugal barrier. 

responsible for the resonances. The bound-state calcu- 
lations used the coupled-states approximation with the 
same basis set as for the collision calculations. The he- 
licity was restricted to K = 1 since this is the only value 
that contributes to the llu —> 111 cross section. The 
resulting bound states are shown in Fig. QTJ Since the 
rotational energy of the bound and quasibound states of 
the Rb-ND3 complex is approximately proportional to 
L(L + 1). it is convenient to plot the energies as a func- 
tion of this quantity. This allows rotational progressions 
to be identified as nearly straight lines, and rotational 
constants can be extracted from the slopes of the lines. 
However, if there are two or more bound states for the 
same L with energies very close together, they can mix 
and avoid one another; an example of this is shown by 
the dotted lines in Fig. [11] 

Most of the peaks in Fig. \lU\ can be identified with se- 
ries of bound states. For example, the cross sections for 
L = 11 shows three strong peaks. The first of these, at 
about 0.5 cm -1 , is a nonrcsonant peak associated with 
the barrier maximum. At higher energies are two res- 
onant peaks, labeled as I, and II in Fig. [10] They are 
associated with characteristic features in the eigenphasc 
sum and can be identified as Feshbach resonances asso- 
ciated with the series of bound states along the dashed 
lines I and II in Fig. [JTJ series II also explains the strong 
peak in the cross section for L = 10 near 0.6 cm -1 . Lines 
I and II both have a large slope, corresponding to a rota- 
tional constant B m 0.04 cm -1 and an effective distance 
between Rb and ND3 around R = 9 ao- It is also not dif- 
ficult to assign the resonances labeled by IV, V, VI and 
VII to the appropriate series of bound states; all these are 
Feshbach resonances that arise from bound states with 
small effective intermolecular distances R < 10 do- For 
all these bound states, inspection of the wavefunction 



near the Van dcr Waals minimum confirms that they are 
dominated by rotationally excited basis functions with 

i>i- 

The levels along the line III in Fig. [11] which give 
strong resonances for L = 8 and 9, are associated with 
quasibound states of significantly smaller rotational con- 
stant, B « 0.008 cm -1 , and an effective intermolecular 
distance around 20 ao- The dominant contributions to 
the wavefunction in this case are the two j = 1, k = 1 
states; the u and I states contribute almost equally and 
with opposite signs, corresponding to a single umbrella 
state of non-inverting ND3. This suggests that these are 
shape resonances, and indeed for L = 8 and 9 they lie 
below the energy of the centrifugal barrier, shown as a 
dashed curve in Fig. QT] For L > 9 the extrapolated 
energy of line III is above the barrier maximum and no 
well-defined shape resonances occur. 

It is important to appreciate that the absolute ener- 
gies of the resonances and bound states close to thresh- 
old are strongly dependent on the potential energy sur- 
face. Thus, the studies reported in this section should 
be treated as a paradigm for understanding the resonant 
behavior of cross sections, rather than as predictions of 
the actual energies at which resonances will appear for 
Rb-ND 3 . 



C. Prospects for sympathetic cooling 

As described in the Introduction, there is great inter- 
est in the possibility of sympathetic cooling, in which 
molecules are cooled by thermal contact with a gas of 
laser-cooled atoms such as Rb. The simplest sympathetic 
cooling experiment would overlap an electrostatic trap 
containing a sample of cold molecules such as ND3 with 
a gas of atoms in a magnetic or magneto-optical trap. 
It is therefore of considerable interest to explore colli- 
sion cross sections for promising candidate atoms and 
molecules in the temperature regime between 1 [iK and 
100 mK. 

An electrostatic trap works only for molecules in low- 
field-seeking states, and for NH3 and ND 3 the low-field- 
sccking states correlate at low field with the upper state 
of the tunneling doublet. It is thus very important to 
know whether collisions with ultracold Rb will cause re- 
laxation to the lower tunneling state, which is a high- 
field-seeking state and cannot be trapped electrostati- 
cally. More specifically, we need to know whether the 
cross section for the relaxation from llu to 111 state is 
sufficiently small compared to the clastic cross section. 

To explore this we have performed close-coupling cal- 
culations for very low energies using the same rotational 
basis set as for the CS calculations described above. In 
Fig. [T2] we show the elastic and total inelastic integral 
cross sections calculated for partial waves J = to 8, 
which gives reasonable convergence for collision energies 
up to 100 mK. The contribution to the total cross sections 
from s-wave scattering ( J = 1 partial wave) is dominant 
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in the fiK regime. The Wigner threshold behavior for 
s-wave scattering is followed for kinetic energies up to 
about 100/xK. Elastic cross sections are larger than in- 
elastic cross sections for collision energies above 100 /iK, 
but never by much more than a factor of 10. For sym- 
pathetic cooling, the commonly stated rule of thumb is 
that the ratio of elastic to inelastic cross sections should 
be at least 100. Thus it appears that the inelastic cross 
sections are too large to allow ND3 molecules in low-ficld- 
sccking states to be cooled to sub-mK temperatures by 
collision with ultracold Rb atoms. 
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FIG. 12: Ultracold limit of the elastic Hit — > llu and inelastic 
llu — > 1 1Z cross sections and the contribution from s-wave 
scattering ( J = 1), from close-coupling calculations. 



There remains the possibility of using sympathetic 
cooling for molecules in high-field-seeking states, which 
can be confined using an alternating current trap [50l.l5ll|. 
Confinement of Rb atoms in such a trap has also been 
demonstrated [52| , but the AC frequencies required are 
quite different in the two cases. It is therefore important 
to know whether sympathetic cooling of ND3 molecules 
in high-field-seeking states might be feasible using mag- 
netically trapped Rb atoms, which are themselves in low- 
field-seeking states that are not the ground state in the 
applied ma gne tic field. In previous work on Rb-OH col- 
lisions [l9l. l53j|. we found large cross sections for low- 
energy collisions that changed the hyperfine state of Rb. 
We initially anticipated [2fj] that this would occur for 
Rb-ND3 collisions as well. However, the more thorough 
analysis of the collision Hamiltonian in section II. B above 
suggests that for molecules in closed-shell singlet states 
the atomic spins are likely to be almost unaffected by 
collisions. We therefore now consider that sympathetic 
cooling of high-field-seeking states of ND3 (or NH3) by 
magnetically trapped Rb (or another laser-cooled atomic 
gas) has a good prospect of success. 



D. Sensitivity to the interaction potential 



Collision calculations at ultralow kinetic energies are 
very sensitive to details of the potential used in the cal- 
culations. Quantitative theoretical predictions of param- 
eters related to scattering cross sections, such as the high- 
est bound states or scattering lengths, are possible only 
for the simplest, lightest systems such as metastable he- 
lium (&41 loa ] or hydrogen- lithium mixtures (56| . The 
calculated Rb-NH3 interaction potential we use in this 
paper suffers from many uncertainties resulting from in- 
completeness of the electronic basis set, the approximate 
treatment of electronic correlation, relativistic effects, 
etc., which limit its accuracy to at best a few percent. 

To explore this sensitivity, we introduce an additional 
scaling factor A into the definition of the interaction 
potential. We varied the scaling factor between 0.90 
and 1.05 and explored the variation of the s-wave elas- 
tic fTnu^iiti and inelastic <tii u _>iii cross sections (using 
close-coupling calculations) for a collision energy of 10 
/iK. The results for a representative slice of the range 
of A studied are shown in Fig. [TjJJ The elastic and in- 
clastic cross sections vary dramatically as bound states 
cross threshold and appear as scattering resonances. In 
the range of A explored, the elastic and inelastic cross 
sections pass through more than ten resonances. There 
are some values of A, such as near A = 0.99, where the 
inelastic cross section <Tn u ->iu is significantly suppressed 
by the presence of a strong resonance, but for most other 
values of the scaling parameter the elastic-to-inelastic ra- 
tio ranges lies between 10 and 0.1. Although it is possi- 
ble that the real potential might produce inelastic cross 
sections low enough to allow sympathetic cooling for low- 
ficld-sccking states, it is quite unlikely 
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FIG. 13: Sensitivity of the elastic llu — > llu and inelastic 
llu — > 111 cross sections for RD-ND3 at a collision energy of 
10 uK to a scaling factor A applied to the interaction potential. 
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IV. CONCLUSIONS 

We have studied collisions between Rb atoms and 
NH3/ND3 molecules, motivated by recent progress (9l.[Tol| 
in collision experiments involving velocity-controlled and 
Stark-decelerated beams of molecules and by interest in 
the possibility of sympathetic cooling. We focused prin- 
cipally on NH3 or ND3 molecules initially in the upper 
component of the inversion doublet for j = 1, k = 1. 
This level correlates with the low-field-seeking state in 
an electric field, which can be slowed by Stark decelera- 
tion and trapped in an electrostatic trap. We considered 
collision energies between and 100 cm -1 . Using the 
coupled-states approximation, we calculated the clastic 
cross sections and state-to-state inelastic cross sections 
from the low-field-seeking j = 1 , k = 1 state to other 
rotation-inversion levels. 

The inelastic cross sections are smaller than expected 
for such a strongly coupled system, but are still only 
about a factor of 10 smaller than the elastic cross section 
over most of the energy range considered. Both the elas- 
tic and inelastic cross sections show dense structure due 
to scattering resonances. The resonances in the clastic 
cross sections are diffuse and rather weak compared to 
the background. The total inelastic cross sections have 
much stronger resonances compared to their background, 
especially at collision energies below about 20 cm -1 and 
one can consider the inelastic scattering as mostly res- 
onant in nature. These resonances are washed out for 
larger collision energies. For energies below 90 cm -1 the 
Hit — > 111 inelastic cross section makes the largest con- 
tribution to the total inelastic cross section. Transitions 



to j = 2 and 3 do not change the total inelastic cross 
section significantly as the new channels become ener- 
getically accessible at higher collision energies. 

We have considered the origin of the scattering reso- 
nances, using calculations of the bound states of Rb-ND3 
near threshold. Since the long-range attraction between 
Rb and ND3 is very strong and the reduced mass is fairly 
large, the resulting centrifugal barriers arc much smaller 
than in lighter systems. Most of the resonances appear- 
ing in the cross sections arise from quasibound states in 
which Rb interacts with rotationally excited NH3, and 
can be classified as Feshbach resonances. 

Finally we studied the ultracold limit of the elastic and 
inelastic cross sections, in order to test whether low-field- 
sceking molecules could be cooled to microKclvin tem- 
peratures by contact with an ultracold gas of Rb atoms. 
Hypcrfinc effects were neglected. We found that inelas- 
tic collisions are strong, and the elastic/inelastic ratio is 
unlikely to be sufficient to achieve sympathetic cooling 
of ND3 molecules in the llu state. However, there is a 
good prospect that sympathetic cooling will be possible 
for high-ficld-sccking states of ND 3 or NH 3 , even if the 
atoms used as a coolant are in low-ficld-sccking states. 
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